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The self-interaction corrected local spin density (SIC -LSD) formalism and the standard GGA treatment of 
the exchange-correlation energy have been applied to study the collapse of the magnetic moment of Fe impu- 
rities in MgO. The system Mgi_vFe,-0 is believed to be the second most abundant mineral in the Earth's lower 
mantle. We confirm the experimentally found increase of the critical pressure upon iron concentration. Our 
calculations using standard GGA for a fixed Fe concentration show that different arrangements of Fe atoms can 
remarkably shift the transition pressure of the high spin (HS) to low spin (LS) transition. This could explain the 
experimentally found broad transition regions. Our results indicate that the HS-LS transition in Mgi_,Fe,0 is 
first order. We find that SIC -LSD fails to predict the divalent Fe configuration as the lowest energy configuration 
and discuss possible reasons for it. 
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I. INTRODUCTION 

The high spin (HS) to low spin (LS) transitions in 
magnesiowiistite, Mgi-vFcj-O, and wiistite, Fei-jO, are of 
great geophysical importance. Magnesiowiistite is believed 
to be the second most abundant mineral in the Earth's 
mantle;ii2 While MgO was thoroughly studied both in the- 
ory and experimen t-^'"^'^'^-^-^'^'"'-"-'^''-^i'"^i'^''^ at high pressures, 
few theoretical studies exist for Mgi-iFcj^O mainly because 
of the difficulty to treat the effects of strong correlation 
for the Fe-3c/ orbitals i'^''^''^ The magnetic phase transition 
could strongly influence the partition coefficient of Fe be- 
tween MgO and MgSiOs perovskite and postperovskitCj^ 
dramatically change radiative conductivities,- and lead to 
a hardening of the materials;2i However, magnesiowiistite 
Mgi_tFe^O and wiistite Fei_^0 have been studied in nu- 
merous experimental papers i^"-^^'^^-^^-^^'^^'^^ The main prob- 
lem has been, that wiistite exists only as a non-stoichiometric 
compound4^i22i22i2L2^ There is a general consensus that its 
rocksalt-type structure contains a fully occupied O^" sublat- 
tice, while Fe exists predominantly in the form of Fe^^, with 
some Fe-'^, as well as vacancies and interstitials, exhibiting a 
short range order The resulting clusters arrange with a long 
range order, which - depending on the density of interstitials 
- can be commensurate or incommensuratci ^^-^^'^^'^^ 

For Fei_vO Jacobsen et al.= reported a transition from a 
cubic to a rhombohedral phase at a pressure between 22.8 GPa 
and 27.7 GPa at room temperature. At low temperatures, on 
the other hand, Struzhkin et al. - found a phase transition from 
a cubic to a rhombohedral structure with a Neel-temperature 
of - 198 K at ambient pressure. At ambient temperature 
this transition is shifted to a pressure of 15 GPa. 

At pressures of 85 to 143 GPa a magnetic high-spin to low- 
spin transition was observed in several experiments! '^i^^-^s 
Some studies find a broad transition region^^ while in others^ 
the transition is first order. Diff'erential stress might have 
smeared out the phase transition over a broad region in some 
studies or small compositional differences within the sample 



could have changed the transition behavior 

Recent developments in high pressure physics allowed to 
investigate the LS-HS transition in magnesiowiistite. Badro 
et al.^ and Lin et al.— could prove what was suggested a 
long time before;^ the complete electronic reaiTangement of 
Fe in Mgi-j^Fe^^O under pressure, while the NaCl-type struc- 
ture remains stable across the phase transition. The transition 
pressures are in reasonable agreement (49 to 75 GPa for x - 
0.15^ 50 to 60 GPa for x = 0.25^). The transition is Hnked 
to a hardening of the materials and a decrease of the molar 
volume (across the phase transition: Vols/^ohs = 0.904, bulk 
modulus: K"^^ = 160.7 GPa, K\;^ = 250 GPa21). 

In another study the HS-LS transition in Mgi-^Fe^^O was 
studied for different iron concentrations (x - 0.2, 0.5 and 
0.8)i^ The transition pressure turns out to depend linearly on 
the Fe concentration and is 40, 60 and 80 GPa, respectively. 
For Feo.gvO the transition pressure is shifted to 90 GPa. At 
high Fe concentration {x - 17%) Badro et al.— in an X-ray 
emission study find that the HS state is stable up to 143 GPa, 
whereas Pasternak et al. '''' conjecture from Mossbauer spec- 
troscopy a strong temperature dependence of the transition 
pressure. They found the transition to be completed only at 
120 GPa (Feo.940) at 450 K. 

In the present paper we study the effect of Fe impuri- 
ties in MgO at geophysically relevant pressures, using self- 
interaction corrected local spin density (SIC-LSD) method 
which allows to treat the localized d electrons of Fe on equal 
footing with the other itinerant electrons. In addition, we use 
GGA approach and the VASP code^'^ to study the influence of 
possible Fe impurity clustering in the MgO supercells on the 
relevant transition pressures. The aim is to realize different 
Fe concentrations in MgO and investigate their effect on the 
experimentally observed properties. 

The outline of the paper is as follows. In Section II, we 
briefly describe the SIC-LSD methodology. Section con- 
centrates on the discussion of the present SIC-LSD and GGA 
results, while Section llVj concludes the paper. 
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II. METHODOLOGY 
A. The self -interaction corrected LSD 

The standard LSD approximation for the exchange- 
correlation energy introduces an unphysical interaction of an 
electron with itself, the so-called self-interaction (SI). It is the 
aim of the self-interaction corrected local spin density approx- 
imation to construct a self-interaction free energy functional 

iiSIC = ^LSD - ^(^a"^' (1) 

a 

where a numbers the occupied orbitals and the self-interaction 
correction (SIC), (J^'^, of the orbital a is 

<5r = f/K] + fix'" (2) 

It is known that the exact exchange-correlation energy 
Exclfia], depending on the spin density «„, for the case of a 
single electron orbital with electron density «„, cancels the 
Hartree energy U[na] identically; 

t/[«J+£xc[n„] =0. (3) 

In the case of approximate exchange-correlation energy func- 
tionals this cancelation can be guaranteed by Eq. [1] 

Varying the above SIC-LSD energy functional with respect 
to the orbital spin densities, with the constraint that the (p^'s 
form a set of orthonormal functions, one gets the SIC-LSD 
generalized eigenvalue equations 

I 0„ > = [Ho^ + y,f ^(r)) I 0„ 2 A„„, I cf>„, >,(4) 

a' 

with Hoa- being the orbital independent LSD Hamiltonian. 
The Lagrangian multipliers Aaa' are used to secure the ful- 
filment of the orthonormality constraint. 

Due to the orbital dependent SIC potential, V^^"^, the SIC 
energy functional is not stationary with respect to infinitesi- 
mal unitary transformations among the orbitals. The so-called 
localization criterion 

< 0/? I ^ - I 0a V(a,/J) (5) 

has to be fulfilled to ensure that the solutions of the SIC-LSD 
equations (|4| are most optimally localized to reach the abso- 
lute minimum of the SIC-LSD functional ([T]). 

The SIC-LSD approach is fully ab initio and introduces no 
adjustable parameters, either for the delocalized (band-like) 
or localized electrons. For extended states the SIC vanishes. 
The SIC-LSD formalism has the advantage that it allows to 
compare different valence states and spin configurations of 
the same atom. The nominal valence, Nyai, in the SIC-LSD 
approach is defined as 

N,„l^Z-Ncore-Ns,C, (6) 



TABLE I: The most relevant spin configurations of Fe, labeled HS 
and LS according to Ref . [s^. Here HS state corresponds to the case 
where SIC is applied to all the majority and one minority electrons. 
For realizing LS state, SIC is applied to the three majority and three 
minority t2g electron states. The HS^ stands for the trivalent Fe con- 
figuration, where all the majority d electrons are treated as localized, 
by applying SIC. The calculations with applying no Sl-corrections 
have been labeled Isd. All possible electronic configurations with 
five or six Sl-corrected orbitals namely trivalent and divalent ions 
have been considered in the present study. However, the ones listed 
below correspond to the lowest enthalpies. Here 1 means that a given 
state has been Sl-corrected, and otherwise. 
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where Z is the atomic number, Ncore is the number of core (and 
semicore) states and Nsic is the number of self-interaction 
corrected states. The ground state valence is the one de- 
fined by the ground state energy. In our application of SIC- 
LSD to (Mg,Fe)0 all possible electron configurations of the 
Fe atom with five and six Sl-corrected d orbitals have been 
considered.— The most energetically relevant ones are given 
in Table H] The energy functional without any SIC states is 
simply equivalent to the standard LSD approximation and la- 
belled Isd in the following. Thus LSD is a local minimum 
of the SIC-LSD functional, corresponding to the case when 
all the electrons are treated as itinerant and described by the 
Bloch wave functions. The choice of the SIC orbital config- 
uration, i. e. the number and symmetry of localized vs. de- 
localized states, can significantly influence the corresponding 
total energy. The electronic configuration giving rise to the 
most negative total energy, with the most optimally localized 
orbitals, defines the absolute energy minimum of the SIC-LSD 
energy functional as well as the corresponding valency. 

The SIC-LSD approach used in this work has been im- 
plemented in the linear muffin-tin orbitals (LMTO) band 
structure method with the atomic sphere approximation 
(ASA)f22il2iii The slightly overlapping ASA spheres approx- 
imate the polyhedral Wigner-Seitz cell and the sum of their 
volumes equals to the volume of the actual unit cell. The 
ASA does not allow for atomic relaxations, as different ionic 
arrangments will lead to different sphere overlaps. In the 
LMTO-ASA approach empty spheres (E) can be introduced 
in order to increase the space filUng and to minimize the over- 
lap. 

In all the calculations performed with the LMTO-ASA 
method we have treated the valence states of the Mg-3s, Mg- 
3p, 0-2s, 0-2p, Fe-3d, Fe-4i, Fe-4p, E-li as the so-called low 
waves, while the Mg-Bc/, 0-3d and E-2p states have been rep- 
resented as intermediate waves 4^ All the electrons in lower 
shells have been put in the core but allowed to relax. The 
separation of valence electrons into the low and intermediate 
waves has been done for the purpose of reducing the size of 
the eigenvalue problem. The secular equation that is fully di- 
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agonalized is constituted by the low waves and provides the 
eigenvalues and eigenvectors. The intermediate waves enter 
the Hamiltonian through their tails, partly retaining their true 
characteristics. 

For both magnesiowiistite and wiistite, we have assumed 
the ideal rock salt structure. The ASA radii of MgO were 
adjusted in order to reproduce the experimental band gaps of 
7.83 eV4^ Best agreement was found with radii of 1.34646 A, 
1.05313 A, 0.751847 A for Mg, O and the interstitial spheres 
E, with the resulting LDA band gap of 4.7 eV. The total vol- 
ume of the spheres was equal to the volume of the cell for all 
the SIC calculations while the ratios of the ASA radii were 
kept constant. The ASA radius of Fe was equal to the one of 
Mg. 

To realize different concentrations of Fe impurities in the 
systems studied, three supercells have been used, namely con- 
sisting of four, eight and 32 formula units of MgO. By sub- 
stituting up to four Mg atoms by Fe atoms one could real- 
ize Fe concentrations (x = -^^) of 3.125%, 12.5%, 25%, 
50% and 100%. For Fe concentrations of 3.125%, 12.5%, 
25%, we have chosen the smallest possible supercell of the 
ideal FCC/NaCl, namely replacing one Mg atom by one Fe 
atom respectively in the supercells consisting of 32, eight and 
four formula units. The concentrations of 50% and 100% 
have been realized using a supercell of eight atoms in to- 
tal and replacing respectively two and four Mg atoms by Fe 
atoms. Consequently, both ferromagnetic (FM) and antiferro- 
magnetic (AFM) orders as well as charge disproportionation 
could be studied. This corresponds to an Fe arrangement with 
maximal distances between the Fe atoms and would provide 
the lowest transition pressures. Regarding the AFM structure, 
the chosen supercells allowed to realize only the so-caUed 
AFl order where the magnetic spins are arranged in parallel 
in the (001) planes, but are anti-aligned between the planes. 
The number of k-points has been chosen inversely propor- 
tional to the size of the supercell (eight, 16, 64 atoms) and 
was 16x16x16, 8x8x8 and 4x4x4, respectively. This 
gives energy differences precise to within 1.7x10"^ eV/atom 
(1.25x10"^ Ryd/atom). 



B. Pressure determination 



The pressures and enthalpies, of relevance to the present 
study, can be determined from the LMTO-ASA total energy, 
E, vs. volume, V, calculations through invoking the equa- 
tion of state (EOS) and fitting to a number of calculated data 
points Specifically, the third-order Birch Murnaghan EOS 
has been used^i^ namely 



16 ^ 

+(x^'^ - if X (6 - 4x^'^)), (7) 
Pi^) - '-^{x^'^-x^'^) 

xjl + ^(B'_4)x(x-/3-l)j, (8) 

X = ^ , (9) 

allowing to compute the enthalpy H at any volume as H{V) - 
E{V) - P{V) V, with P denoting the pressure. In the above, Vq 
stands for the equilibrium volume, Bq for the bulk modulus 
and B' for its first pressure derivative. 

The calculation of enthalpy differences, AH(P), can be ac- 
complished through the numerical inversion of the pressure in 
Eq. |8] This allows to compare enthalpies of two spin configu- 
rations labelled s and s' at a given pressure: 

AH'-"'^ = H\V') - W\V'') . (10) 

To determine the theoretical equilibrium volumes and total 
energies for all the studied scenarios, the lattice parameters 
have been sampled in the range from - 10% up to + 16% (rel- 
ative to the zero-pressure value) in steps of 2%. The resulting 
volume increase or decrease of -33.1 % to 40.7% typically cor- 
responds to a pressure range of -20 GPa to 285 GPai^ 

III. RESULTS AND DISCUSSION 
A. (Mg,Fe)0 in the SIC-LSD framework 

In the application of SIC-LSD to Fe-doped MgO, we have 
studied 29 different localization-delocalization scenarios for 
three different supercells, various Fe concentrations, and sym- 
metries of localized Fe d states. Among all the cases we have 
identified four energetically relevant configurations, namely 
the three SIC configurations listed in Table I] as well as the 
LSD solution. This finding has been independent of the actual 
Fe concentration. The HS configuration is the high-spin diva- 
lent scenario, where all majority Fe d states and one minority 
Fe t2g state are treated as localized through invoking SIC, giv- 
ing rise to an insulating state. The HS^*, trivalent, configura- 
tion corresponds to the scenario where only the majority Fe d 
states are corrected for the self-interaction, and thus localized, 
constituting a half-metallic state. The LS state corresponds 
to the case where both majority and minority Fe tog states are 
locaUzed by SIC. This naturally leads to a non-magnetic and 
insulating configuration. 

From the enthalpy calculations we find, in disagreement 
with experimental evidence, that the HS^ and LSD solutions 
are energetically most favorable, respectively at the low and 
high pressure regions. This can be seen in Fig. |l(a)[ where 
the enthalpy differences, with respect to the HS"* configura- 
tion, are plotted as a function of pressure for all the relevant 
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scenarios and Fe concentration of 3.125%. The trivalent HS"' 
state constitutes the ground state solution for all the pressures 
up to about 140 GPa. However, at very low pressures the diva- 
lent HS state lies very close in energy to the HS^ state and the 
enthalpy difference, at GPa and 3.125% Fe concentration, is 
^^HS - ^^HS' =0.05 eV. Note that at this Fe concentration the 
total energies of the HS and HS"' configurations, evaluated at 
their respective theoretical volumes, are fully degenerate. 

The nominal valence of the HS"* state, being 3+, is at odds 
with the experimental findings, where Fe appears clearly as 
Fe-^ in (Mg,Fe)0 (see e.g. Ref. l20h or with a small ratio of 
additional Fe^^ in FexO (see e.g. Ref. [32l) . Note, however, 
that the nominal valence defined in SIC-LSD not always cor- 
responds to the chemical valence. Although the nominal va- 
lences of the trivalent and divalent states differ by one (namely 
3+ vs. 2+), in terms of a simple charge counting or charge 
disproportionation, the two ions, corresponding respectively 
to the HS^^ and HS configurations, differ only by up to 0.1 
electron. Similar observation was made by Szotek et al.— for 
magnetite. Despite the latter, it is still surprising to find the 
trivalent, instead of the divalent, Fe-state as the lowest en- 
ergy configuration, in particular for Fe^^O (with x=l). This is 
in variance to the SIC-LSD results obtained for all the other 
transition metal monoxides, namely MnO, CoO and Ni04i^ 
It is possible that the failure of SIC-LSD to find the divalent 
ground state in FeO has its origin in the fact that LSD sub- 
stantially overestimates the exchange splitting for the systems 
in question. Consequently, the energy gained on localizing an 
additional electron to create a divalent Fe ion is not sufficient 
to overcome that exchange splitting. However, since FeO does 
not occur in a stoichiometric form, but as Fei_ vO, and exper- 
iments indicate existence of both divalent and trivalent ions, 
one would need to perform more realistic calculations to es- 
tablish the importance of the off-stoichiometry and valence 
fluctuations for obtaining the divalent ground state in FeO. 

At high pressures, it is the experimentally indicated-- LS 
configuration that is relevant, with its nominal valence of 2+ 
in agreement with the experimentally found Fe-^ in wustite 
and magnesiowiistite. The LS configuration competes with 
the LSD description of the Fe-3d states, as seen through the 
corresponding enthalpy difference which at GPa and 3. 125% 
Fe concentration is //kd - //ls =0.7 eV (see Fig.lTJ, in favour 
of the LS state, while the LSD becomes energetically more fa- 
vorable for pressures in the excess of 140 GPa. At such high 
pressures the gain in band formation energy on delocalization 
of all the Fe d electrons clearly wins with the localization en- 
ergy of those electrons. Since the magnitude of SIC is strongly 
dependent on the orbitals, one needs to be guided by energet- 
ics in defining the ground state energy and configuration. 

Introducing several Fe impurities into a supercell 
(Mg2Fe204 and Fe404) opens a possibility of realizing 
both parallel (labelled [f], ferromagnetic) and antiparallel 
(labelled [a], antiferromagneticsAFl) arrangements of the 
spin magnetic moments on those Fe atoms. In Fig. |l(b)| the 
enthalpy diff'erences are plotted for the above mentioned con- 
figurations, and both parallel and antiparallel arrangements 
of spins on the Fe atoms in the relevant supercell with 50% 
Fe concentration. For the HS configuration the [f] and [a] 
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FIG. 1: Enthalpy differences at different Fe concentration x. |(a)| 
X = 3.125%. Tfie HS^ configuration appears to be energetically most 
favorable up to a pressure of 145 GPa, where a phase transition is 
predicted to the state where all the Fe-3d delocalize (Isd). The com- 
monly assumed HS configuration is less energetically favourable, 
and at the pressure of 90 GPa the electrons completely rearrange to 
the non-polarized configuration LS, although the localization of the 
electrons is not lost. \(b)\ x = 50%. Two Fe atoms were introduced 
into a cell containing 8 atoms in total (two Mg, two Fe, four O). Par- 
allel and antiparallel alignments of the spin moments is denoted in 
brackets [f] and [a], respectively. More than 10 different arrange- 
ments with asymmetric charge configurations on the Fe atoms were 
also tried out, but resulted in high enthalpies. They are omitted in the 
plot for more clear readability. 



arrangments return the same energies while for HS the [f] 
arrangement appears always lower in energy {AE -0.5 eV 
for Mgo.5Feo.5O and Afi =1.5 eV FeO*"^) than the antiparaflel 
alignment. 

Finally, we consider the dependence of the HS-LS transi- 
tion pressure on the concentration of the Fe impurities in a 
supercell. As can be seen in Fig. |2l our study predicts an 
increase of this transition pressure with rising Fe concentra- 
tion. However, the absolute transition pressures for isostruc- 
tural phase transitions, as calculated in this SIC-LSD study, 
are higher than the experimental values42 Two possible sce- 
narios can operate at any Fe concentration: a possible transi- 
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FIG. 2: Mgi_vFejO: Experimental and theoretical transition pres- 
sures. Squares (■) - SIC-LSD calculations, scenario HS'-lsd (elec- 
tronic configuration HS^^ at low pressures, no Sl-corrections at 
high pressure); bullets (•) - SIC-LSD calculations scenario HS- 
LS (electronic configuration HS at low pressures, configuration 
LS at high pressure); lozenges (♦) - GGA-PAW calculations (low 
concentrations); stars (★) - LDA-l-U calculations tsoll ; triangles 
(a) - Mossbauer spectroscopy, (36ll ; arrow (J) - Mossbauer spec- 
troscopy (3^ ; circle (O) - X-ray diffraction (32I1 . 



tion from HS^ to Isd or from HS to LS. The second one occurs 
at considerably lower pressures than the first one, and is the 
relevant one for comparison with experiments, although lead- 
ing to a phase diagram in rather poor agreement with experi- 
ment. The latter implies the existence of Fe-^, and a low pres- 
sure structure with one minority Fe V^g state and five majority 
d stateSf2S, One might envisage that the HS state could possibly 
become favourable if in the SIC-LSD calculations the ionic re- 
laxations had been applied. It is obvious that ASA error can 
substantially affect total energies, although one would like to 
hope that these errors would be less severe when considering 
only the energy differences between different configurations. 
Also, across the phase transition Fe is known to reduce its 
ionic radius and therefore the ionic relaxations could signifi- 
cantly reduce the enthalpies of the LS structures which in turn 
could further decrease the transition pressurei^ 



B. Density of states from SIC-LSD 

In Fig. [3] we present the densities of states (DOS) calcu- 
lated for all the energetically relevant scenarios and the Fe im- 
purity concentration of 12.5%. As can be seen in the figure, at 
this intermediate concentration the localized Fe-3ii states form 
narrow bands. At higher concentrations of Fe impurities these 
bands extend over a considerable energy range, e.g. 3-5 eV at 
25% of Fe. For the HS and HS^ configurations, the occupied 
Fe "id band states lie below the predominantly O 2/7 valence 
band. As mentioned earlier, the SI coiTections for the HS and 
HS^ configurations differ by one localized electron more in 
the former, which is reflected in the presented DOS (see Figs. 
|3(a)| and |3(b)| i. The difference is that the first isolated, minor- 
ity Fe t2^ band, lying just below the valence band, is moved 
above the valence band, coinciding with the Fermi level (see 
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FIG. 3: |(a)| Spin-resolved density of states (DOS) for HS state as cal- 
culated from SIC-LSD at a pressure of 2 GPa and Fe concentration 
of 12.5%. |(b)| Spin-resolved density of states for the HS^ state as 
calculated from SIC-LSD at a pressure of 2 GPa and Fe concentra- 
tion of 12.5%. |(c)| Spin-resolved density of states for the LS state 
calculated at a pressure of 124 GPa and Fe concentration of 12.5%. 
|(d)| Spin-resolved density of states calculated with the standard LSD 
approximation at 102 GPa and Fe concentration of 12.5%. In all the 
cases the Fe DOS states is in red, while the total DOS is in black and 
the O 2/7 DOS in green. 



Fig. [3(b)l ). While for the HS configuration SIC-LSD dehv- 
ers, in agreement with with experiment and other theoretical 
considerations ,^ ''^^'^-^ an insulator of charge transfer character, 
the HS^ configuration gives rise to a half-metal, with a large 
band gap in the majority band, and a metallic behavior in the 
spin-down channel, with a partially occupied minority Fe d 
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band at the Fermi level. 

At high pressures the LS and LSD configurations are of rel- 
evance. The DOS of the LS is shown in Fig. |3(c)| It is charac- 
terized by hybridized Fe-3ii electrons with the predominantly 
2p valence band, and a large band gap to the Fe d conduc- 
tion band. The Isd DOS shows a substantially reduced band 
gap, which is between the occupied and unoccupied Fe im- 
purity d bands. So, while the SIC-LSD LS state is a charge 
transfer insulator, the LSD gives rise to a Mott-Hubbard insu- 
lator (see Fig. |3(d)| i. 



C. GGA-PAW calculations, Fe clustering 

So far we have concentrated mostly on the correlated na- 
ture of Fe d electrons in describing electronic properties of 
doped MgO within SIC-LSD approach. Here we consider 
an importance of a possible Fe impurity clustering. To study 
this, we have performed a number of independent calculations 
for low Fe concentrations in Mgi_ rFe^O, using the standard 
GGA approximation-'* and the Vienna Ab-initio Simulation 
Package (VASP)^^*^. The latter provides the standard projec- 
tor augmented wave potentials for Mg, Fe and O with core 
configurations 2p''3s- (Mg), 3p^3d^4s' (Fe) and 2s-2p'* (O)^ 
In the actual calculations the electronic optimization has been 
run with a self-consistency threshold of 10"' eV, ionic op- 
timization has been done with a threshold of 10"^ eV. The 
4 X 4 X 4-Monkhorst-Pack scheme has been applied for the 
k-point sampling, producing four irreducible k-pointsi^ The 
energy cut-off of the plane wave basis set has been equal to 
700 eV. In order to speed up the evaluation of the non-local 
part of the potentials real space projections have been used. 
A second order Methfessel-Paxton smearing has been used 
(cr=0.04 eV), introducing an error in electronic entropy, T ^ei, 
of at most 0.111 me V/atom^^ The accuracy of the energy dif- 
ferences between the HS and the LS state has been converged 
to 1.3 X 10"^ eV/atom. Note, however, that within the GGA 
implementation all the electrons are treated on equal footing 
as delocalized, with the HS state being described as a ferro- 
magnetic state, whilst the LS state as a non-magnetic state. 
Thus the two are not directly comparable to the HS and LS 
states calculated within SIC-LSD, as the concept of divalent 
and trivalent ions cannot be clearly defined within GGA. 

To investigate the effects of Fe clusters, supercells contain- 
ing 64 atoms have been set up, including one and four Fe im- 
purities (xpe - 3.125 %, xpe =12.5 %). In the latter case, 
30 distinct arrangements of Fe atoms exist3 We have used 
the most extreme cases, where the Fe atoms form a tetrago- 
nal cluster with an Fe-Fe distance of ^ 2.69 A and where 
they avoid each other at a maximal distance of ^ 5.37 A (at 
100 GPa). The resulting three clusters (one Fe atom, four Fe 
atoms at minimum distance, four Fe atoms at maximum dis- 
tance) have been explored performing spin-polarized and non- 
spin-polarized GGA calculations (in the former case, a start- 
ing magnetic moment of 4 //^ has been assigned to each Fe 
atom). In the spin-polarized case the Vosko-Wilk-Nusair in- 
terpolation for the correlation part of the exchange-correlation 
functional has been usedj^. The lowest energy structure has 



been found by comparing the enthalpies of all the different 
structures. 

The GGA-PAW calculations in the 64 atomic supercell 
find HS-LS transition pressures of 26 GPa (x - 3.125 %), 

30.3 GPa {x = 12.5 %, diluted configuration of Fe), and 

63.4 GPa (x = 12.5 %, clustered configuration of Fe), as 
calculated from the intersection of the enthalpy curves of 
the ferromagnetic and non-magnetic results. The lowest en- 
thalpy structures always give a total magnetization of either 
fi = O.O/Yfi (non-magnetic) or ju - A-.QjUB (ferromagnetic). Ar- 
rangements of Fe atoms with an intermediate spin magnetic 
moment have been found in some cases (e.g. fj. - 3/iB in 
the case of Fe clusters between 90 GPa and 185 GPa). How- 
ever they always turn out to be unstable because energetically 
lower lying cluster configurations exist. 

The transition pressure turns out to increase with the Fe 
concentration. The two studied concentrations {x - 3.125% 
and X = 12.5%) allow to determine the linear dependence of 
the transition pressure as 



Per = 24.6 + 0.45 X . (11) 

The relatively small difference to the experimental value of 
Per = 28+0.63 X can be explained, if it is assumed that the sta- 
tistical distribution of the Fe atoms in the experimental sample 
will not be perfectly diluted, but in some cases result in having 
Fe atoms at smaller distances (which dramatically raises the 
critical pressure). However, this can also explain the broad- 
ness of the experimentally found transition pressure: depend- 
ing on the local Fe configuration the spin magnetic moments 
collapse at different pressures. 

Figure |4]indicates that at low pressures (where only the HS 
arrangments are of interest) the Fe atoms prefer to form clus- 
ters, while at high pressures (where only the LS arrangments 
are of interest) the Fe atoms prefer a diluted configuration. 
Upon pressure release at the pressure of about 60 GPa the LS 
diluted arrangement and the HS clustered arrangement coin- 
cide. The different Fe arrangments in the high pressure and 
low pressure regions might hamper this phase transition. Ex- 
perimental HS-LS transition pressures are predicted to depend 
on the pressure the sample has been equilibrated at (the equi- 
libration could be enhanced through heat, which accelerates 
ionic diffusion). If the equilibration is done at high pressure, 
upon pressure release the LS configuration would be stable up 
to relatively low pressures (^30 GPa). A sample equilibrated 
at low pressure on the other hand would preserve its magneti- 
zation up to relatively high pressures (^60 GPa). 



IV. CONCLUSIONS 

Our GGA results indicate a first order transition for all the 
compositions studied, both for Mgi_vFexO and Fei^^O, in 
good agreement with the notion in Ref. [33 for Fei_ vO. The 
electronic rearrangement of the Fe atoms leads to a discon- 
tinuity of the first derivative of the enthalpy, i.e. the volume 
(see Fig.|4]i. This affects various physical properties of the HS 
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HS, diluted 



HS, clustered 




LS, clustered 



MgO+4 Fe, 64 atoms cell comparision 



P [GPa] 



FIG. 4: GGA-PAW calculations, 64 atomic supercell including four Fe atoms: Comparison of the enthalpies of the different Fe and spin 
arrangements to the diluted LS setting. In the clustered configuration Fe atoms are at the largest distance possible. 



and LS phases, such as the zero pressure bulk modulus and 
the ground state volume. 

In experiment the HS-LS transition appears to be smeared 
out over a large pressure range,is2i21^ which in the case of 
Fei_vO has led to the conjecture that the transition might be 
second order^ We can not support this unless this is a tem- 
perature induced smearing which our T = calculations do 
not consider. The results for various Fe arrangements, keep- 
ing the Fe concentration fixed, indicate that the local arrange- 
ment of Fe strongly influences the spin transition pressure 
of the Fe atoms. The statistical distribution of the Fe atoms 
and the short range defect order can therefore smear out the 
transition pressure over a large pressure range. As the HS- 
LS transition is isosymmetric, according to Landau theory it 
can only be first-order or (above the critical temperature) fully 
continuous 

We can confirm the experimentally found trend that the HS- 
LS transition pressure increases with increasing Fe concentra- 
tion. This might be the result of Fe atoms which — if situated 
closely enough — simultaneously optimize their electronic 
configurations. These magnetic Fe clusters can not easily be 
demagnetized by pressure. The increase of the Fe concen- 
tration leads to a statistical increase of short Fe-Fe distances, 
which through the described mechanism could in turn lead to 
an increase of the critical pressure. 

It has been proposed that the HS-LS transition in magne- 
siowiistite would lead to a discontinuity in the Earth's lower 
mantle;^ Although the transition appears to be first order a 
discontinuity seems to be improbable: first the sensitivity of 
the transition pressure to the (local) Fe concentration could 
smear out the transition over a wide region in depth. Sec- 
ond, temperature could further increase this effect and shift 
the transition to even higher pressures, barring LS Fe^^ from 
the Earth's interiorJ^i^ 

The transition pressures predicted by SIC-LSD appear to be 
higher which might be a result of an "overcorrection" of the 
energy of the localized states. It has been argued in Ref. 
that a weighted self-interaction correction would be desirable. 
The high transition pressures calculated in this work might 



also be due to the fact that the ASA does not allow relaxations 
of the atomic spheres, which might be significant in the LS 
phase, where the ionic radius of Fe is small. At high pressures 
the LSD treatment can lead to low energy structures. This 
could justify the treatment of Mgi-vFCfO at those pressures 
using standard approximations for the exchange-correlation 
energy such as the GGA. Apart from this it could also indicate 
a second phase transition which involves a full delocalization 
of the Fe-3ii states. Furthermore all the transitions considered 
in this paper involve five or six localized electrons in the high 
volume phase and no localized electrons in the low volume 
phase. Consideration of valence fluctuations or 'intermediate' 
localization involving four, three, two and one localized states 
could plausibly lead to a reduction of the transition pressure. 
The SIC-LSD formalism, implemented in the multiple scat- 
tering theory, allows to consider static valence and spin fluc- 
tuations at finite temperatures, by invoking CPA (coherent po- 
tential approximation) and DLM (disordered local moments) 
approaches, which was successfully applied to studying phase 
diagramsi^^ This would also be useful for the further theo- 
retical investigation of the present geophysically relevant sys- 
tems. For the nearest future, it may be instructive to first apply 
SIC-LSD in a full potential version in order to allow for ionic 
relaxations to study the influence of distortions on the ground 
state properties of these systems. 



Acknowledgements 

The authors are indebted to A.R. Oganov who contributed 
considerably to the concept of this work. DJA thanks A.R. 
Oganov and the ETH Zurich Research fund for their sup- 
port of this work (Grant No. TH-27033). Supercomputers 
were provided by the Swiss National Supercomputing Centre 
(CSCS) and ETH Zurich. WMT and ZS would like to ac- 
knowledge useful discussions with Drs. Axel Svane and Leon 
Petit. 



8 



* Electronic address: donat.adams @ cea.frl 

t CEA, DAM, DIP, F 91297 Arpajon, France 

1 S. Ono and A. R. Oganov, Earth Planet. Sci. Lett. 236, 914 (2005). 

^ D. L. Anderson, Theory of the Earth (Blackwell Scientific Publi- 
cations, Oxford, 1989). 

^ A. R. Oganov, Ph.D. thesis, University of London (2002). 

" A. R. Oganov and P L Dorogokupets, Phys. Rev. B 67, 224110 
(2003). 

^ T. S. Dufi'y, R. J. Hemley, and H. kwang Mao, Phys. Rev. Lett. 74, 
1371 (1995). 

J. E. Jaffe, J. A. Snyder, Z. Lin, and A. C. Hess, Phys. Rev. B 62, 
1660 (2000). 

^ H. Zhang and M. S. T. Bukowinski, Phys. Rev. B 44, 2495 (1991). 

^ K. J. Chang and M. L. Cohen, Phys. Rev. B 30, 4774 (1984). 

' B. B. Karki, R. M. Wentzcovitch, S. de Gironcoli, and S. Baroni, 

Phys. Rev. B 61, 8793 (2000). 
1° N. D. Drummond and G. Ackland, Phys. Rev. B 65, 184104 
(2002). 

" G. Kalpana, B. Palanivel, and M. Rajagopalan, Phys. Rev. B 52, 4 
(1995). 

'2 M. Causa, R. Dovesi, C. Pisani, and C. Roetti, Phys. Rev. B 33, 
1308 (1986). 

" A. Strachan, T. gagin, and W. A. Goddard, Phys. Rev. B 60, 15084 
(1999). 

A. Dewaele, G. Fiquet, D. Andrault, and D. Hausermann, J. Geo- 
phys. Res. 105, 2869 (2000). 

S. Speziale, C.-S. Zha, T. S. Duffy, R. J. Hemley, and H. Mao, J. 
Geophys. Res. 106, 515 (2001). 

M. Mehl, R. E. Cohen, and H. Krakauer, J. Geophys. Res. 93, 
8009 (1988). 

K. Persson, A. Bengtson, G. Ceder, and D. Morgan, Geophys. 
Res. Lett. 33 (2006). 

T. Tsuchiya, R. M. Wentzcovitch, C. R. S. da Silva, and S. de 

Gironcoh, Phys. Rev. Lett. 96, 198501 (2006). 
" Z. Fang, I. V. Solovyev, H. Sawada, and K. Terakura, Phys. Rev. 

B 59, 762 (1999). 
2" J. Badro, G. Fiquet, F Guyot, J.-P Rueff, V. V. Struzhkin, 

G. Vankao, and G. Monaco, Science 300, 789 (2003). 
^' A. F. Goncharov, V. V. Struzhkin, and S. D. Jacobsen, Science 

312, 1205 (2006). 
22 E. Gaffney and D. L. Anderson, J. Geophys. Res. 78, 7005 (1973). 
2^ Y. Fei, H. kwang Mao, J. Shu, and J. Hu, Phys. Chem. Minerals 

18, 416 (1992). 

2" J.-F Lin, V. V. Struzhkin, S. D. Jacobsen, M. Y. Hu, P Chow, 
J. Kung, H. Liu, H. kwang Mao, and R. J. Hemley, Nature 436, 
377 (2005). 

2^ T. Welberry and A. G. Christy, Phys. Chem. Minerals 24, 24 
(1997). 

26 T. Welberry and A. G. Christy, J. Solid State Chem. 117, 398 
(1995). 

2^ C. Carel and J. Gavarri, Phys. Chem. Minerals 26, 78 (1998). 
2** J.-R. Gavarri and C. Carel, Rev. Chim. Minerale 18, 608 (1981). 
2' Y. Ding, J. Xu, C. T. Prewitt, R. J. Hemley, H. kwang Mao, and 

J. A. Cowan, Applied Phys. Lett. 86 (2005). 
^° Y. Ding, H. Liu, J. Xu, C. T. Prewitt, R. J. Hemley, and H. kwang 

Mao, Applied Phys. Lett. 87, 041912 (2005). 
^' Y. Ding, H. Liu, M. Somayazulu, Y. Meng, J. Xu, C. T. Prewitt, 

R. J. Hemley, and H. kwang Mao, Phys. Rev. B 72, 174109 (2005). 
32 S. Ono, Y. Ohishi, and T. Kikegawa, J. Phys. Cond. Matt. 19, 

036205 (2007). 

'3 S. D. Jacobsen, J.-F. Lin, R. J. Angel, G. Shen, V. B. Prakapenka, 



D. Przemyslaw, H.-K. Mao, and R. J. Hemley, J. Synchroton Ra- 
diat. 12, 577 (2005), international workshop on Structure Deter- 
mination by Single-Cristal X-ray Diffraction (SXD) at Megabar 
Pressures, Chicago , USA (13/1 1/2004). 

V. V. Struzhkin, H. kwang Mao, J. Hu, M. Schwoerer-Bohning, 
J. Shu, R. J. Hemley, W. Sturhahn, M. Y. Hu, E. E. Alp, R Eng, 
et al., Phys. Rev. Lett. 87, 255501 (2001). 

M. P Pasternak, R. D. Taylor, R. Jeanloz, X. Li, J. H. Nguyen, 

and C. A. McCammon, Phys. Rev. Lett. 79, 5046 (1997). 
36 S. Speziale, A. Milner, V. E. Lee, S. M. Clark, M. P Pasternak, 

and R. Jeanloz, Proc. Natl. Acad. Sci. 102, 17918 (2005). 
" J. Badro, V. V. Struzhkin, J. Shu, R. J. Hemley, and H. kwang 

Mao, Phys. Rev. Lett. 83, 4101 (1999). 

G. Kresse and J. Furthmuller, Phys. Rev. B 54, 11169 (1996). 
3' W. M. Temmerman, A. Svane, Z. Szotek, , and H. Winter, in 

Electronic Density Functional Theory: Recent Progress and New 

Directions, edited by J. F. Dobson, G. Vignale, and M. P. Das 

(Plenum, New York, 1998), pp. 327-347. 
*° O. K. Andersen and O. Jepsen, Phys. Rev. Lett. 53, 2571 (1984). 

O. K. Andersen, Phys. Rev. B 12, 3060 (1975). 
■•2 W. R. L. Lambrecht and O. K. Andersen, Phys. Rev. B 34, 2439 

(1986). 

*^ R. C. Whited and W. C. Walker, Phys. Rev. Lett. 22, 1428 (1969). 
J.-P. Poirier, Introduction to the Physics of the Earth's Interior 
(Cambridge University Press, Cambridge, UK, 2000). 

■•^ F Birch, Phys. Rev. 71, 809 (1947). 

■•6 Z. Szotek, W. M. Temmerman, A. Svane, L. Petit, G. M. Stocks, 
and H. Winter, Phys. Rev. B 68, 054415 (2003). 
D. Kasinathan, J. Kunes, K. Koepernik, C. V. Diaconu, R. L. 
Martin, I. D. Prodan, G. E. Scuseria, N. Spaldin, L. Petit, T. C. 
Schulthess, et al., Phys. Rev. B 74, 195110 (2006). 
M. Dane, M. Liiders, A. Ernst, D. Kodderitzsch, W. M. Temmer- 
man, Z. Szotek, and W. Hergert, journal of Physics: Condensed 
Matter 21, 045604 (2009). 

A. Svane, P. Strange, W. M. Temmerman, and Z. Szotek, Phys. 
Stat. Sol. 223, 105 (2001). 

''° T. Tsuchiya, R. M. Wentzcovitch, C. R. S. da Silva, S. de Giron- 
coli, and J. Tsuchiya, Phys. Stat. Sol. B 243, 2111 (2006). 

^' F. Parmigiani and L. Sangaletti, journal of Electron Spectroscopy 
and Related Phenomena 98-99, 287 (1999). 

52 G. A. Sawatzky and J. W. Allen, Phys. Rev. Lett. 53, 2339 (1984). 

" J. Zaanen, G. A. Sawatzky, and J. W. Allen, Phys. Rev. Lett. 55, 
418 (1985). 

5'* J. P. Perdew, K. Burke, and M. Emzerhof, Phys. Rev. Lett. 77, 
3865 (1996). 

" G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 

56 H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 (1976). 

" M. Methfessel and A. T. Paxton, Phys. Rev. B 40, 3616 (1989). 

^'^ S. Vosko, L. Wilk, , and M. Nusair, Can. J. Phys. 58, 1200 (1980). 

5' A. G. Christy, Acta Cryst. B51, 753 (1995). 

6" A. D. Bruce and R. Cowley, Structural phase transitions (Taylor 

& Francis, London, 1981). 
6' W. Sturhahn, J. M. Jackson, and J.-F. Lin, Geophys. Res. Lett. 32 

(2005). 

62 J. P. Perdew, A. Ruzsinszky, J. Tao, V. N. Staroverov, G. E. Scuse- 
ria, and G. I. Csonka, J. Chem. Phys. 123, 062201 (2005). 

63 M. Liiders, A. Ernst, M. Dane, Z. Szotek, A. Svane, 
D. Kodderitzsch, W. Hergert, B. L. Gyorffy, and W. M. Temmer- 
man, Phys. Rev. B 71, 205109 (2005). 

6* I. D. Hughes, M. Dane, A. Ernst, W. Hergert, M. Liiders, J. Poul- 
ter, J. B. Staunton, A. Svane, Z. Szotek, and W. M. Temmerman, 



9 



Nature 446, 650 (2007). 

The nominal valence A'vai = Z - A'core - ^sic does not always 
correspond to the chemical valence. For example, in the case of 
the conventional LSD scheme the valence of Fe would be A'vai = 
Z-A'co,e-A'sic=26-18-0 = 8. 

We have used a standard procedure such as provided with the EX- 
CITING code at http://exciting.sourceforge.net 
The enthalpy differences have been fitted to the second order poly- 
nomials. However, the error resulting from these fits and the nu- 
merical inversion of Eq. [8] have given rise to total errors of the 
order 10"** - 10"'^ eV/atom. 

AE is given for one supercell containing 2 Mg, 2 Fe and 4 O in 
the first case and 4 Fe and 4 O in the second. 



In our calculations the Fe-AS A radii were fixed to those of Mg at 
ambient conditions. According to Anderson— they are compara- 
ble: rMg = 0.72 A, rpc.HS = 0.77 A. However, in the LS phase the 
radius of Fe dramatically reduces to 0.61 A. Due to the LMTO- 
ASA approach, the reduction of the ionic radii could not be taken 
into account. 

Note that this is a combinatorial problem. The clusters can be 
characterized by the distance between the first Fe impurity to the 
following ones. Five different distances can be found in the 64 
atomic cell. They can host 12, three, 12, three and one Fe atom 
respectively. This gives one, five, 14, 30, 53 distinct scenarios of 
one, two, three, four and five impurities, respectively. 



